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ABSTRACT 

During the inspiral and merger of a binary black hole, gravitational radiation is emitted anisotrop- 
ically due to asymmetries in the merger configuration. This anisotropic radiation leads to a gravita- 
tional wave kick, or recoil velocity, as large as ^ 4000km s^^. We investigate the effect gravitational 
recoil has on the retention of intermediate mass black holes (IMBH) within Galactic globular clusters. 
Assuming that our current understanding of IMBH-formation is correct and yields an IMBH-seed in 
every globular cluster, we find a significant problem retaining low mass IMBHs ( ^ 1000 M©) in the 
typical merger-rich globular cluster environment. Given a uniform black hole spin distribution and 
orientation and a Kroupa IMF, we find that at most 3% of the globular clusters can retain an IMBH 
larger than IOOOMq today. For a population of black holes that better approximates mass loss from 
winds and supernovae, we find that 16% of globulars can retain an IMBH larger than 1000 Mq. Our 
calculations show that if there are black holes of mass A/ > 60 Mq in a cluster, repeated IMBH-BH 
encounters will eventually eject a, M = 1000 Mq IMBH with greater than 30% probability. As a 
consequence, a large population of rogue black holes may exist in our Milky Way halo. We discuss 
the dynamical implications of this subpopulation, and its possible connection to ultraluminous X-ray 
sources (ULXs). 

Subject headings: black hole physics — galaxies: nuclei — gravitation — gravitational waves — 
relativity 



1. INTRODUCTION 

Ample observational evidence exists for two 
types of black holes: stellar mass (BH), 
with lOM© < m < 10^ Mq, and supermass ive (SMBH), 
with m > 10^ Mq (e.g.. lKormendv fc Ricl^tond (|1995,1). 
Although the existence of a third black hole class is 
still under debate, there are observational hints for 
intermediate mass black holes (IMBHs) as well, with 
masse s lO'^ Mc, < m < Mrr. fe.g.. iGebhardt et ^ 
(l2005l): iGerssen et alJ (l2Q0l: iFilippenko fc Hd 



clusters may require more than one single mechanism. 
We briefly introduce these mechanisms below. 

Recent simulations of core collapse in stellar clus- 
ters have shown that stellar collisions could induce 
rapid stellar runaway gro wth of an 0(10'^ Mq) IMBH 
in a. s little as 3 Myr (iPortegies Zwart fc McMillanI 



iHeggie et all (l2006bD: iTrenti et all (l2007dlal) 
' 2006D : lUlvestad et all (|2007| ): c.f. iBaumgard 



2003li 



Trenti 



et al 



2003| ) for an alternative view.) These IMBHs may 
form within dense star clusters and may be the 
best explanation for ultraluminous X-ra y sources 
(ULXs) in young star forming regions (iFabb land 
19891 [Roberts fc Warwicli [2000I: iPtak fc ColbertI 120041 : 



Fabbiano fc Whitd I2006D and in nearby extragalactic 



star clusters. 

If the observational case for IMBHs is still unclear, the 
exact formation mechanism is perhaps less constrained. 
Within a globular cluster, there are currently three plau- 
sible IMBH formation theories: stellar runaway; compact 
object merger s; and growth from Population (Pop) HI 
remnants (cf. Ivan der Marell (|2004f) for a review.) Each 
of these mechanisms starts with very different initial con- 
ditions, and the first two require the high stellar densi- 
ties of a globular or stellar cluster environment. Thus, 
growing IMBHs over the ensemble of observed globular 
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12002"; 'Giirkan et al.' '2006'; Freitag_eialJ 120061); see also 
(Portcgies Zwart ct al. 2004). However, this runaway 
growth process requires number densities higher than 
O(10^/pc'^) and few (if any) globular clusters in the 
Milky Way are this dense today. Of course, globu- 
lar clus ters are thought to be more dens e at formation 
(jTrenti et al. 2007b; He ggie et all l2006aD . but calcula- 
tions still suggest that the runaway process may only 
occur in ab out 20% of the Milky Way globular clus- 
ter system (jBaumgardt et al.l I2005D . By using less ex- 
treme central densities and including the effect of pri- 
mordial binaries, simulations have shown that instead of 
a single IMBH, stellar runaways g enerate two or mor e 
black holes with masses ~ 500Mq (iGiirkan et al.l[200il . 
When these large black holes merge they form a sin- 
gle IMBH and a strong gravitational wave signal for the 
Laser Interferometer Space Antenna, or LISA, a planned 
gravita tional wave observatory set to launch in the next 
decade ()Fregeau et al.ll200l IGiirkan et aTll2006( ). 

While runaway star collisions seem to be the preferred 
IMBH formation channel within dense systems, in sys- 
tems where the relaxation time is long, stellar mass BH 
collisions may be a more likely growth mechanism. Here 
the initial seed IMBH forms naturally from a heavier 
than average bl ack hole, perhaps a supernova remnant 
with ~ 250 Mq (|Wisefc Abell 120051 ). Then, stellar mass 
BHs are captured by the initial seed, and the resulting 
binary hardens through interactions with other stellar- 
mass objects. Although these interactions risk ejecting 
the BH supply before significant growth occurs, one way 
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to funnel more BHs to the growing IMBH is through 4- 
body colHsions (lO'Learv et alJ [20061 : iMiller fc Hamiltonl 
[200a iGiiltekin et al.ll200 

Significantly higher mass Pop III stars may yield 
IMBHs directly. If a Pop III protostar has a mass of 
IO^Mq, it will be gravitationally unstable and will col- 
lapse to an IMBH before it even enters the main se- 
quen ce (jBaumgarte fc Shapirolll999l : IShibata fc Shapirol 
120021 ). The probability of an IMBH forming from such 
a massive Pop III star certainly depends on the initial 
stellar mass function (IMF), as well as highly uncertain 
details of zero-metallicity stellar evolution. However, 
there has been some suggestions that the IMF in the 
early universe is quite top-heavy ([Schneider et al.l [20021 : 
[Abel et aD '20021. and that stellar mass loss is negligible 
([Fryer et al . 2001; Heg cr et a l. 2003). These suggestions 
imply that BHs in a proto globular cluster environment 
are more massive tha n those formed in more recent times 
([van der Marelll2004t ). On the other hand. Pop HI star 
formation is thought to take place in dark matter over- 
densities at z 12 — 20 (Madau fc R ccs 2001) and only 
the most massive globular clusters are thought to be em- 
bedded in a significant enough dark matter overdensity 
to allow for Pop HI formation. 

However the IMBH forms, for the fi rst ^ 0.5 Gyr after 
formation it is dynamically active (Spitzer ,1987). Due 
to mass segregation, the initial environment around an 
IMBH in a globular cluster is especially rich in BHs 
([Fregeau et all l2002t l. While many of these BHs are 
quickly ejected by few body interactions with the IMBH, 
enough remain to subject the IMBH to 10 s — 100 s merg- 
ers with BHs in the primordial globular cluster systeni 
(iPorteeies Zw art fc McMillanll2000HO'Learv et al.ll2006l : 
IGiiltekin e t al. 2003). 

In light of general relativistic black hole merger simu- 
lations, surviving this IMBH-BH merger epoch may be 
difficult. Recent advances in numerical relativity have 
at last pinned down the dynamics of black hole merg- 
ers - simulating the coalescence, merger, and ringdown 
of equal-mass circular non-spinning binary black holes in 
full g eneral relativi ty (jP rctorius 2005; Campanell i et al.) 
[20051: [Baker et al.[ [200l . In the past year, more as- 
trophysically relevant numerical simulations, including 
spins and unequal masses, have been published by sev- 
eral groups. 

One of the most exciting results of general relativity 
for structure formation is that binary black hole sys- 
tems strongly radiate linear momentum in the form of 
gravitational waves during the plunge phase of the inspi- 
ral. This radiation directly results from an asymmetry 
in the orbital configuration and can generically yield a 
gravitational wave "kick" velocity as fast as 4000 km 
([Gonzalez et~al[[2007bt ICampanem et al.l[2007b[ ). Even 
typical kick velocities (~ 200 km s~^) are interestingly 
large when compared to the escape velocity of an av- 
erage globular cluster (^ 50km s~^) ((WebbinkI Il985[ : 
[Favata et all [20041 : [MCTritt et all [2004 ). Hence, regard- 
less of how an IMBH forms, the biggest challenge may 
be to determine how a globular cluster retains them in 
the face of a repeated onslaught of gravitational wave 
kicks from mergers with other black holes. Of course if 
IMBHs are formed through mergers of stellar-mass black 
holes, the formation mechanism itself needs to be able to 



account for these large kicks as well. 

In this paper we explore IMBH retention within young 
globular clusters after collisions with BHs. We calcu- 
late the retention probability under a variety of assump- 
tions for the initial IMBH seed mass, the BH mass dis- 
tribution, and the initial spin distributions. In addition, 
for the mass range encompassed by the IMBH formation 
channels listed above (stellar runaway, stellar-mass BH 
collisions. Pop HI stars), we simulate the IMBH survival 
probability after merging with the available BHs over a 
collision timescale. Given the known current structure 
of the Milky Way globular cluster system, we can then 
estimate the fraction of globular clusters that may have 
allowed a particular IMBH formation channel. Combin- 
ing these two key estimates, we can determine the max- 
imum expected number of Milky Way globular clusters 
that could have retained their IMBHs. We outline the 
approach in Sec. 2, present the results in Sec. 3 and dis- 
cuss the caveats, implications, and future directions in 
Sec. 4. 

2. ASSIGNING KICKS 

Gravitational recoil estimates of binary black 
hole mergers have been addressed using both 
semi-analytic metho ds in the u nequal-mass, 
non-s pi nning case ([Fitchettl [19831 : [Favata et al.l 
'2004- 'Damour fc Gopakumad [20061 : [Blanchet et"an 
2005; Sopuerta et al. 200^ and numerical meth- 
ods in mo re general, sp inn ing, un equal mass 
scenarios ([Herrmann et al. | 20071 : [Baker et al 



[Gonzalez et aL [2006: [ Herrmann et al 
'Koppitz et all [2007t [G onzalez et al.' '2007a[: 
Campanelli et al. 2 007a[ [d: [Tidw fc M arronctti 2007,). 
For non-spinning binaries, the most comprehensive nu - 
merical study of kicks is that of [Gonzalez et ahl ([2006f) . 
which was fo und in agreement w ith the semi-analytic 
estimates of ([Sopuerta et al.ir2006( ). where a maximum 
kick velocity of ~ 175 ± 10 km s~^ was obtained for 
mass ratio q = M1/M2 0.36 ± 0.02). 

The gravita tional recoil is expected to increase with in - 
creasing spin ([Redmount fc R,eesl[l98l lWhitbeckll2(¥)i) . 
and this behavior has been confirmed by several numeri- 
cal relativity groups. Numerical simulations have demon- 
strated that the radiative linear momentum loss pre- 
dicte d by Post-Newtonian studies ([Thornel[1980t [Kiddeil 
[T99I can be used to fit numerical results, yielding a gen- 
eralized formula for the recoil velocity as a function of the 
individual black hole's spin, initial orientation, phase at 
merger a nd mass ratio. W e here a dopt the parameter- 
ized fit of [Campanelli et all ([2007a[ ) addi ng the expected 
(1_+ e) contribution for eccentric orbits ([Sopuerta et al.l 
[20071 ) to yield the following formula: 

Vkick = (1 + e) [x {vjn + v± cos^) + yvj_ sin^ -I- z vn] , 

(1) 

where 



A 



(1 - q) 



l + B- 



(1 + 9) 



H- 



(2) 



(3) 



and 



if cos (9 - Oc 



{a^ - qai) , (4) 
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where the fitting constants are A = 1.2 X lO'^ km 
B = -0.93, H = (7.3 ± 0.3) x 10^ km s"\ and K = 
(6.0 ± 0.1) X lO'* km s^^, while the subscripts 1 and 2 
refer to the first and second BH respectively. The unit 
vectors (i,y) are orthogonal to each other and span the 
initial orbital plane, while _L and || stands for perpendicu- 
lar and parallel to the orbital angular momentum. There 
are 4 fitting parameters: the mass ratio q = M2/M1; the 
reduced spin parameter ai = Si/Mf, where Si is the 
spin angular momentum of BH z; and the eccentricity 
e. In addition, there are 3 angles to specify the orienta- 
tion of the merger: 0, the angle between the "in-plane" 
component oi 5' = {M1+M2) {S\/M2 - S\/Mi) and the 
infall direction at merger; Go, the angle between 5^ and 
the initial direction of motion; and ^, the angle between 
the unequal mass and spin contribution to the recoil in 
the orbital plane. The merger phase within the orbital 
plane may also play a role, but it is not included explic- 
itly in this fit. The recoil velocities as given in Eq. ([4]) 
are plotted in Fig. [T] as a function of mass ratio and spin 
parameter. 

In order to compute the probability that an IMBH re- 
mains in its globular cluster, our calculation shall pro- 
ceed as follows. First we begin by assuming that an 
IMBH has formed within a globular cluster with a par- 
ticular initial mass, Mimbh, which we shall vary. Sec- 
ond, we further assume a certain mass distribution for 
the BHs in the vicinity of the IMBH, which we shall also 
vary. Third, we subject this IMBH to a number of merg- 
ers expected within a proto globular cluster environment 
that we shall describe below. Finally, we determine the 
probability that the kick velocity for the IMBH has re- 
mained below the canonical globular cluster escape ve- 
locity (50km s^^) during the entire chain of mergers. 

Even in the absence of an IMBH, BHs eject themselves 
from globular clusters via standard few-body interactions 
on a timescale of ~ 1 Gyr after the onset of mass segrega- 
tion (iKulkarni et al .' 1993; Sigurdsson & Hcrnquist 1993; 
iPortegies Zwart fc McMillan 2000 : O' Leary et al. 2006). 
Therefore, due to such Newtonian few-body interactions, 
the supply of BHs is eventually depleted. With an IMBH, 
however, this process speeds up impressively, as ejec- 
tions by interactions with an IMBH become the domi- 
nant source of stellar-mass BH ejections (Giiltek in et al.l 
I2OO6D . As in most few-body interactions, the ejection of 
one object tightens the orbit of a remaining bound pair, 
in this case an IMBH-BH binary ~ and after several sub- 
sequent ejections, the IMBH-BH binary merges. Soon 
after all the BHs have been evacuated, the short epoch 
of IMBH-BH mergers ends. 

Within this theoretical framework, it is possible to con- 
struct a fiducial number of mergers for a proto globular 
cluster. This number can be written as 



''^cjcctions/mcrgcr 



(5) 



Gultekin et al. 2006 predict iVmergo ~ 25 per IMBH, 
and we adopt this for the fiducial number of mergers 
that the IMBH encounters. Although we do vary this 
parameter in figure [21 the dependence of the retention 
probability on the number of mergers is relatively minor, 
since the IMBH grows in mass over each merger and the 
kick velocity increases with increasing mass ratio. 



In order to assign a kick velocity to each of the mergers, 
we choose the orientation, spin, mass, and eccentricity. 
We outline the assumptions made for each distribution 
below. Let us first discuss the issue of the initial spin 
orientation. Hydrodynamic interactions between a gas 
disk and a black hole binary are believed to align the spin 
directions to the angular momentum axis of the binary 
orbital plane in many active galaxies (jBogdanovic et al] 
l2007f ). However, the environment of a globular cluster is 
not particularly gas-rich, so there is no ab initio reason to 
expect the black hole spins to be aligned. We therefore 
assume an isotropic distribution of orientation angles for 
each encounter. 

Let us now discuss the choice of spin magnitude. Most 
theories predict a non-zero spi n for a black hole p ro- 
duced via stellar runaway (,Rees fc Volonter"i| | 2007f) or 
from a supernovae remnant ( Frver et al.l l200l[ r If an 
IMBH started with zero spin, a merger is likely to spin 
up the remnant through transfer of o rbital to spin angu- 
lar momentum (|Gammie et al.|[2004D . However, a Kerr 
black hole can spindown when magnetic field lines thread 
through the ergosphere to m agnetically brake the system 
(|Blandford fc Znaieklll977[ ). provided there is a gaseous 
disk around the remnant. Taking all these considera- 
tion into account, we shall explore three cases: (1) the 
spin magnitude is selected from a uniform initial spin 
distribution (fiducial case); (2) the initial spin magni- 
tude of the IMBH seed is set to 0.998Mj2^bjj; and (3) 
the spin is initially set to zero. We assume the spin of 
the secondary BH to be randomly selected from a distri- 
bution of [0,0.998] M^^^, where Msec is the mass of the 
secondary BH. Since these stellar-mass BHs originate as 
a sup ernova remnant, t hough, the spin is likely to be 
high (jHeger et al.ll2003l ). Therefore, since the kick ve- 
locity increases with increasing spin, this shall lead to a 
conservative survival probability for each merger tree. 

The survival probability shall be studied as a func- 
tion of the initial mass of the seed IMBH. We choose a 
range of 10 — 3OOOA/0 to encompass the plausible IMBH 
formation channels and relevant masses involved. For 
the three formation channels discussed in § [1] the seed 
masses are likely to be t he following: ^ lOOOMf?) fo r 
a single stellar runaway (jPortegies Zwart et al.l 1200^ : 
~ IOOMq for the stellar mass BH collision channel, where 
the seed IMBH is produced by a massive supernova rem- 
nant or a small stellar runaway (He ger e t al. 2003); and 
^ 200 — 400 M(7) for Pop III remnants within a dense glob- 
ular cluster (jWise fc Abel 120051 : iMadau fc Rees 20011 

One of the biggest uncertainties is determining a 
proper distribution for the secondary BH masses. As fig- 
ure [2] shows, the retention probability depends strongly 
on the mass ratio between the IMBH and BH. Theo- 
retical black hole mass distributions from solar metal- 
licity fi eld populations tend to peak strongly around 
IOM0 (jFrver fc Kalogeral [2001h . There are, however, 
several strong competing effects in a primordial globular 
cluster that can change this distribution (e.g., low metal- 
licity stellar evolution, high binary fraction, stellar colli- 
sions, mass segregation, and natal kicks). For our fiducial 
experiment, we assume that the secondary masses are se- 
lected fr om a Kroupa IMF with an upper mass cutoff of 
I2OM0 (|Kroupall200l[ l. Recall that this yields an aver- 
age stellar mass of about IM0, much smaller than the 
suspected average stellar mass from a zero-metallicity 



Fig. 1. — Left: The gravitational wave kick velocity from a non-spinning black hole merger as a function of mass ratio and eccentricity. 
The white c ontour mar ks 50 km s~^, the escape velocity of a typical globular cluster. Overlaid are data from numerical relativity merger 
simulations l IGonzalez et al. 2006.^ (black dots on zero-eccentricity line.) Note that most of this plane has not been confirmed by direct 
numerical experiment, but has been extrapolated from semi-analytic studies (Sopuerta ct al. 2007). In this plane, only mergers with the 
largest mass ratios will remain in the globular cluster, which argues for a very high-mass initial seed. Right: The gravitational wave kick 
velocity from spinning black hole mergers on a circular orbit as a function of spin magnitude and mass ratio. Here, the spins of the black 
hole are anti-aligned to one another and are perpendicular to the orbital plane of the binary black hole before merger. Note that other 
spin orientations yield even higher kicks. As in the previ ous figur e, the 50 k m s~^ contour is o verplotted, as wel l as se lected results from 
full numerical re l ativit y simulations: squares are from Koppitz et al.i (2 0071 ) : circles are from [Herrmann et all (1200711 . and star is from 
IBriigmann et al.l l l200l) . This fi gure suggests that globular clusters can retain IMBHs only when they merge with other BHs with low spin 
and high mass ratio. 



environment. We further assume that each star above 
IOMq evolves directly to a BH with no mass loss. This 
simplified treatment gives us an average BH mass of 
~ 20Mq. Clearly, mass loss would decrease the average 
BH mass even in these low metallicity primordial globu- 
lars. Therefore, we used a more sophisticated black hole 
mass function that approx imates the results of model 
CI of iBelczvnski et al.l (2006). This model includes the 
effect of mass loss from supcrnovac and winds for a pop- 
ulation of stars with metallicity Z — 0.0001, though it 
has fewe r massive remnan ts from binary mergers (cf Fig- 
ure 8 of IBelczvnski et al] ((2006)) than is expected for a 
primordial globular cluster. 

Figure O demonstrates the effect on the kick velocity 
distribution between these two black hole mass functions 
for a IOOOMq IMBH. Since the distribution of secondary 
masses is so uncertain, we demonstrated the effect of 
varying the mass ratio in Figure [H For a near-solar 
metallicity stellar cluster, the entire p opulation of black 
holes may be less than M ^ 20 Mq (jFrver fc Kaloger"al 
12001 . 

Finally, we must select the eccentricity distribution. 
If these IMBH-BH mergers were 2-body processes, we 
might expect the eccentricity of the orbit to be very 
nearly circular right before the black holes merge as 
gravitational r adiation grinds away the orbital angu- 
lar momentum (jPeters fc Mathewill963f) . However, few 
body encounters are much more common within a glob- 
ular cl uster because the inte raction cross-section is much 
larger (jHeggie fc Hutll2003f ). Therefore, many BHs are 
shepherded into mergers with an I MBH through ex- 
chang e with lower mass black holes (jMiller fc Hamilton! 
I2OO2D , and the r esulting eccentricity can be quite 
large (jCiiltekin et al. 2006). In fact, simulations have 
shown that rare interactions can yield mergers with 
e > 0.999, and such a highly eccentric orbit can be- 
come even m ore eccentric through gravitational radia - 
tion emission (jPeters fc Mathewslll963l:lKenneficklll998D . 



Therefore, though rare, highly eccentric binary black hole 
mergers can take place in astrophysically relevant sys- 
tems. To assign eccentricities t o each merger^ we use 
the simulation results of Giiltekin et al.l (|2006[ ). which 
empirically characterizes the eccentricity distribution as 
a function of the mass ratio of the encounter. Note, 
though, that equation ([1]) is really only valid in the small 
eccentricity regime and, thus, it may be true that the 
kicks can be even higher than those studied here for such 
nearly radial orbits. 

3. IMBH SURVIVAL AND OCCUPATION FRACTION 

In order to determine the probability that an IMBH 
survives the short merger epoch after formation, our sim- 
ulations consist of 10^ Monte Carlo realizations of an N- 
step merger chain. The merger chains are tailored to 
mimic the initial conditions and encounters predicted by 
current IMBH formation theories within proto-globular 
and stellar cluster environments as described in the pre- 
vious section. 

During a merger, gravitational waves radiate not only 
linear momentum, but also angular momentum and en- 
ergy or mass. Fully relativistic numerical simulations 
suggest that ^ 25% of the angular momentum (defined 
at the innermost sta ble cir cular orbit) can be radiated 
during the merger (iPretor ius 2005; Campanclli et al.) 
I2005t iBaker et al.l l2006t iHerrmann et al.. .2007. ). In 
addition, the mass of the merger product is only 
~_95%the mass of the two progenito r black holes 
(IPretorius! 2005: Campanein et all 12005!: iBaker et all 
12006!: IH errmann et al.l l2007l ). Hence, after each step 
within a particular merger tree, we adjust the total spin 
and mass of the merger remnant to account for these 
losses. 

Figure Ej demonstrates the retention percentage after 
25 mergers with black holes selected from a Kroupa IMF 
as a function of the initial IMBH mass. This figure in- 
dicates that retaining an IMBH of less than IOOOM0 oc- 
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Fig. 2. — Left: Retention probability as a function of mass ratio of the interacting black holes. We assume an initial IMBH mass of 
A^iMBH = 500 Mq , an escape velocity of 50 km , and a merger chain of 25 interactions. Mergers are chosen with random orientations and 
spin magnitudes. Middle: Retention probability as a function of the number of BH mergers. We assume Mbcc = 20 Mq , AfiMBH = 500 Af , 
and random orientations and spin magnitudes. Right: Retention probability as a function of the escape velocity of the host structure. We 
assume Msec = 20 Mq, Mimbh = 500 Mq, and random orientations and spin magnitudes. 
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Fig. 3. — Fraction of encounters as a function of kick velocity 
for two black hole mass functions. We assume the IMBH has a 
mass of IOOOA/q and that the spin magnitudes and orientations 
are distributed uniformly. 

curs less than 33% of the time for the given distribu- 
tion of black hole masses. While the large number of 
the lower-mass black holes dominate the total number 
of mergers, the rare mergers with massive stellar-mass 
black holes dominate the ejections - and this trend only 
strengthens as the IMBH mass increases. For example, 
figure [5] shows that for a 1000 Mq seed, nearly all ejec- 
tions come from black holes with mass M > 30 Mq 
and that most come from those with mass M ~ 70 M0 . 
Naturally, this implies that if > 3OM0 black holes are ex- 
tremely rare in primordial globu lar clusters, th e retention 
fraction increases dramatically ( Fryer fc Kalog era 2001), 
making gravitational recoil ineffective in ejecting mas- 
sive IMBHs. Figure [6] shows that it is easier to retain an 
IMBH of less than lOOOM© with the shallower Belczynski 
black hole mass function, with only 30% ejected. 

Recall that each formation channel is expected to take 
place in a particular proto-globular cluster environment 
(e.g., single runaways ought to form in very dense sys- 
tems.) Since predictions of pre-core collapse clusters in- 



FlG. 4. — A histogram of the black hole masses that eject 
a IOOOMq IMBH. In this experiment, we assumed a Kroupa 
IMF, and a uniform spin orientations and magnitude. Ejecting 
a IOOOA/q black hole is possible with masses as low as 20AfQ , but 
the overwhelming majority of ejections result from black holes with 
masses > 50Mq. 

dicate that the densities are ^ 10 times higher than to- 
day, we can us e current obs ervations of the central lumi- 
nosity density (|Harrislll996f ) to infer the initial conditions 
of the globular cluster. Then, we can estimate on a case- 
by-case basis what IMBH formation channel might have 
been appropriate for each globular cluster. If we assume 
that every globular cluster forms an IMBH with a mass 
that is consistent with its presumed formation channel 
(see previous section), we can estimate the number of 
surviving IMBHs within the Milky Way globular cluster 
system. We find that ^ 24% of the Milky Way globular 
clusters can host single runaways, while 25% may host 
multiple runaways. If, thereafter, the IMBH must sur- 
vive 25 mergers selected from a Kroupa IMF, our study 
indicates that no more than 5 globular clusters have re- 
tained an IMBH, while our Belczynski black hole mass 
function yields < 25 IMBH-embedded globulars. 

Such conclusions can be split into the number that sur- 
vive per formation channel. Doing so, we learn that the 
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Fig. 5. — Percentage of black hole retention as a function of ini- 
tial black hole mass within the Milky Way globular cluster system. 
A black hole is defined as 'retained' if it survives 25 collisions with 
a BH selected from a Kroupa IMF with random spin orientations. 
The three lines on this figure represent different assumptions for 
the spins of the black hole. 
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Fig. 6. — Percentage of black hole retention as a function of 
initial black hole mass within the Milky Way globular cluster sys- 
tem. In this experiment, we assumed a black ho le mass function 
described by model CI of IBelczvnski et al.l 1)20061 ). and a uniform 
spin orientations and magnitude. 

only way to retain IMBHs during the short BH merger 
phase is with one rather massive initial seed, such as 
those produced in a stellar runaway; any process that re- 
lies of the formation of ^ 500Mq black holes will be inef- 
fective, as these are ejected from a globular cluster before 
forming an IMBH. One way to lower the retention mass 
threshold would be to increase the globular cluster escape 
velocity. We explore this effect in figure [U where we ob- 
serve that the escape velocity must reach ~ 80 km s~^ 
to attain complete retention of 500Mq black holes. 

A BH mass function that is tipped more strongly to 
lower masses would also increase the retention rate - note 



that a Kroupa IMF with no mass loss yields an average 
BH mass of ~ 20Mq, while allowing for mass loss will 
make high mass ratio mergers less likely. If there are no 
black holes > 20AIq in clusters, then any IMBH seed 
larger than M ^ 600 Mq will remain. This critical mass 
ratio is comparable to that necessary to avoid ejection 
from three-body dynamic a l kicks alone (iGiiltekin et al.l 
[200l lO'Learv et all [20061 : [Giiltekin etalj 120061) . Si^e 
runaway stellar mergers preferentially include the clus- 
ter's heaviest stars, which are the precursors to the heav- 
iest black holes in the cluster, the largest black holes used 
in our calculations may not be present in the cluster. 

4. CONCLUSION 

Our studies indicate that it is a challenge to retain 
IMBHs during the barrage of BH mergers expected in an 
early globular cluster. We find no scenario that guaran- 
tees 100% retention for IMBHs up to 3000 Mq when the 
interacting BHs are selected from a reasonable distribu- 
tion of spins, orientations, and mass ratios. However, the 
more massive the initial IMBH seed, the better the reten- 
tion probability. This may indicate that any IMBH ob- 
served in globular clusters today would most likely have 
originated from an early stellar runaway channel. 

The results obtained here might be used to test 
some IMBH forma tion mechanisms. For example, 
iFregeau et all (|2006l ) and lGiirkan et al.l l|2006l ) find that 
with a binary fraction above 10% a single cluster can 
host multiple runaways, leading to multiple IMBHs. The 
heaviest two such IMBHs will find their way to the cen- 
ter of the cluster and merge very quickly ('~ 1 Myr). 
Because such IMBHs are nearly equal in mass and likely 
have high spins, they will escape the cluster upon merger. 
This leaves the cluster without a seed unless a third run- 
away of sufficient mass occurs. The problem is that the 
third runaway is often much less massive and therefore 
difficult to retain, as we have shown. 

In fact, since lower mass IMBHs so easily escape glob- 
ular clusters, if globular cluster observations find that 
they are not rare, it may be possible to constrain the 
IMBH merger history, as well. For example, if many low 
mass IMBHs exist within clusters, we may rule out low 
mass ratio mergers - this would imply that there is no 
high mass tail in the BH IMF. Alternatively, the lower 
the spin, the better the retention as figure [5] shows; if low 
mass IMBHs are found in large numbers within clusters, 
and if BHs are found with > 20Mq , we may have to ex- 
plore ways in which the black holes spin down and align 
within a gas-poor globular. 

Making some very simple assumptions for the primor- 
dial globular cluster environment, such as the BH IMF 
and central density structure, we have estimated that less 
than 5 globular clusters retain their IMBHs within the 
Milky Way - even if every one hosted an initial IMBH 
seed. Naturally, there are many uncertainties folded into 
this estimate, such as the shape of the primordial glob- 
ular cluster IMF, the degree of mass loss in low mctal- 
licity systems and its effect on the BH IMF, and the de- 
tailed role that few-body/BH interactions play in shaping 
early globular cluster structure. As more work is done 
on these areas, we can revise our estimates using more 
generic results. For example, if the BH IMF were in- 
stead narrowly peaked around IOMq (jFrver fc Kaloger"al 
l2001h and if the seed black hole were a massive super- 
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nova remnant, our calculations indicate that about 25% 
of the Milky Way globular clusters could retain black 
holes as small as ~ 400Mo. An IMBH of a few hun- 
dred Mq may not leave an electromagnetically observ- 
able impact on the surrounding globular cluster, as the 
dynamical effects on the surroundi ng stars may also b e 
produced by a high binary fraction (|Trenti et al.ll2007"bl ) . 
However, when stellar mass compact objects merge with 
these smaller mass IMBHs, they will produce a strong 
gravitational wav e signals that s hould be detectable with 
Advanced LIGO (|Mandel et al.l[2b07i) . 

With so many small black holes having been ejected 
from their host globular clusters, we speculate that 
100 rogue black holes are swarming about in the Milky 
Way halo with masses from ~ 100 — IOOOMq, and with 
velocities mostly on the order of a few hundred km s~^. 
The number of rogues could be in the thousands if, as has 
been suggested, the current globular cluster population 
is a small fraction o f the total number originally created 
(|Aguilar et al.lll988f ). 

Although we have focused on IMBHs in the Galactic 
globular cluster population, the same processes may oc- 
cur in other galaxies. Extragalactic ULXs, which may be 
powered by ~ 10^ M© IMBHs, are frequen tly found near, 
but not in you ng stellar clusters (e.g., Fabbian o et al.l 
120011: iLiu fc"Bregman,2005.1 . Note, though, that the stel- 
lar clusters associated with ULXs are not always the 
dense stellar systems re quired by the I MBH formation 
models considered here (|Liu et al.ll2007f) . If such IMBHs 
did form within the nearby clusters, they may be ejected 
from gravitational wave kicks coming from mergers with 
stellar-mass black holes, especially as they would merge 
with the most massive black holes first. While this would 
explain their separation from the cluster center, it would 
not explain the fact that ULXs are accreting sources; it 
is not clear how an IMBH would pick up a companion 
on its way out of the cluster and is unlikely to retain a 
stellar companion close enough to overfill its Roche lobe. 
An IMBH with a stellar companion, however, may be 
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ejected from the host cluster through few-body Newto- 
nian dynamica l kicks that harden the binary until it be- 
gins accreting ( iGiihekin et al.ll2004!:t OT.earv et al.ll2006t 
Giih ekin et all 2006; Blecha et al. 200^^ 

Even if the ejected IMBH is not accreting g ULX, 
electromagnetic observations may still detect rogue black 
holes. For instance, if the IMBH were to carry a few 
massive stars along as it is ejected, our results indicate 
a kinematically fast subpopulation of massive stars near 
globular clusters. The ejected black holes may leave a 
temporary imprint on the globular cluster as well. Since 
they are ejected from the system impulsively, it is likely 
that the globular cluster core would temporarily expand. 
Direct simulations remain to be done to determine how 
the globular cluster responds to the ejection of an IMBH. 

The consequences of these large recoil velocities may 
also affect SMBH assembly. The most likely candidates 
for SMBH seeds are ~ Pop III stellar remnants at 

redsh ifts z > 12 - 20 flleger et al.l 120031: IVolonteri et al. 
'2003; Isla m et al.l [20031: iWise fc Abell 120051: iMicic et al. 
i200 7). These relic seeds are predicted to form at the cen- 
ters of low mass dark matter halos (~ 4 x IO^Mq). As 
dark matter halos hierarchically merge to assemble the 
galaxy, the seed black holes sink to the center through 
dynamical friction and eventually merge. With kick ve- 
locities in the range of ^ 10^ — lO^km s~^, it may also 
be difficult to retain seed SMBHs in high redshift low 
mass dark matter halos. We plan to explore black hole 
retention and possible kick suppression mechanisms at 
the low mass end of the halo mass function using high 
resolution cosmological N-body simulations in our next 
paper. 
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